Spin torques due to diffusive spin current in magnetic texture 
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We present a microscopic theory of spin torque due to diffusive spin currents induced by spin 
accumulation. The obtained expression is a natural extension of the existing one due to 'local' spin 
currents associated with ordinary electric currents, and is the reciprocal of the spin motive force 
which induces charge accumulation as studied recently [J. Shibata and If. Kohno, Phys. Rev. B84, 
184408 (2011)]. The result is applied to a domain wall motion in a nonlocal spin injection system, 
and the torque and force due to diffusive spin current are evaluated. 
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I. INTRODUCTION 

The control of the magnetization by using flow of the 
spin of conduction electrons continues to be one of major 
research topics in the field of spintronics since the pio- 
neering works of Slonczewski^ and Berger^ on multilayer 
structures. Fundamental physics behind these phenom- 
ena is understood as due to spin torque, which arises from 
transfer of the spin angular momentum between con- 
duction electron and magnetization. The idea was first 
applied by Berger to the current-induced domain wall 
motion'^"^, where the electron spin can follow the local 
magnetization when traversing the domain wall. Since 
the theoretical details of the current- induced spin torques 
were exposed by phenomcnological^'^ and microscopic* 
analyses, it is well-known that there are reactive torque 
(also known as spin-transfer torque)^"^^ and dissipative 
torque (also called /3 term)^''"'"'^^"^^ in the presence of spin 
relaxation of conduction electrons. They are expressed 
by the first and the second terms of 



current were carried out^*'^^. The spin torque due 
to diffusive spin current has been discussed theoreti- 
cally in ferromagnet /nonmagnet/ferromagnet / nanopil- 
lar systems^'^'^"'^ and the contribution of the diffusive spin 
current has been estimated quantitatively^^'^''. However, 
a general expression of the spin torque induced by diffu- 
sive spin currents in magnetic texture, such as magnetic 
domain walls, has not been derived microscopically. 

As a reciprocal to the current-induced spin torque, a 
spin motive force is known to be generated by the dy- 
namics of magnetization texture^^"^^ , and actually has 
been detected experimentally^^. In this context, two of 
the present authors studied microscopically the spin and 
charge transport induced by magnetization dynamics^^. 



They obtained a 'local' charge current, j 
driven by a spin motive force field 
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{n ■ {diU X h) + /3ri • diu} , 



(2) 



where the first and the second terms correspond to the 
respective terms in Eq. (1). Moreover, they obtained a 



^ |(is •V)n + /3nx (j^ ■\/)nj, (1) diffusive charge current, arising from charge im- 



respectively. Here n is a unit vector representing the 
direction of the localized spin (magnetization), Js^'^'' = 
((T-|- — o'i)E is a spin-polarized current (spin current) 
driven by a local electric field E, with (cr^) being 
spin-resolved conductivity for majority-spin (minority- 
spin) electrons, and /3 is a constant originating from spin 
relaxation. 

Though it is well recognized that the spin torque 
is induced by the spin polarized current, the case of 
diffusive spin current has not been addressed in pre- 
vious studies. A diffusive spin current is induced 
from spin accumulation and is distinguished from spin- 
polarized current induced by local electric field. Ex- 
perimentally, the local spin-polarized current can be 
separated off from the diffusive spin current by the 
nonlocal spin injection technique^''. By using this 
method, magnetization switching by pure spin current 
and depinning of domain walls assisted by diffusive spin 



balance induced by inhomogeneous Eg. This diffusive 
current; in turn, implies the existence of spin torques in- 
duced by the diffusive spin current as the reciprocal, and 
it is expected that this torque takes the form 



A{(j(D).v)n + /?nx(ip).V)n}, (3) 



where jii^' is the diffusive spin current. 

The purpose of this paper is to clarify microscopically 
the effects of diffusive spin current on the spin torque. 
It is shown that the general expression of the current- 
induced spin torque is given by the sum of Eq. (1) and 
Eq. (3). The result is applied to a domain wall motion 
in a spin-valve system with nonlocal spin injector, and 
the force and torque acting on the domain wall due to 
diffusive spin current are evaluated. 

The paper is organized as follows. After describing 
the model in Sec. II, we calculate in Sec. HI the spin 
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torque in response to an external electromagnetic field. 
In Sec. IV, we apply the obtained result to a domain wall 
motion in a nonlocal spin valve structure and evaluate the 
diffusive spin current and the torque and force acting on 
the domain wall. Conclusion is given in Sec. V. Details 
of the calculation are given in Appendices. 



in the rotated frame is then given by L = L^i — He- a, 
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j dr yl°A>1"a, (10) 



II. MODEL 

We consider a ferromagnetic conductor described by 
the s-d model which consists of conducting s-electrons 
and localized c?-spins. The Lagrangian for s-clcctrons is 
given by 

f f d \ 
Lni = I dr cU ih— + —V^ + ep - Mmp j c - H^d, 

(4) 
(5) 



= —M J dr n{r) ■ {c'crc)^ 



where c^(.t) = (cj(a;), c|(a;)) is the electron creation op- 
erator at X = {t, r), ep is the Fermi energy, M is the s-d 
exchange coupling constant, n(r) is a static but spatially- 
varying unit vector representing the direction of d-spin, 
which is treated as a classical spin, and cr is a vector 
of Pauli matrices. We assume that s-electrons are inter- 
acting with both non-magnetic and magnetic impurities, 
whose potential is modeled by 



^.J2Sr^Sir-Rr), (6) 



where u (ug) is the strength of normal (magnetic) im- 
purities, which introduce momentum relaxation (spin- 
relaxation) processes, and Ri (Rj) is the position of nor- 
mal (magnetic) impurities. 

To treat magnetization texture, we introduce a local 
spin frame for conduction electrons and take the d-spin 
direction n(r) as the spin quantization axis.^^''^^''^'^ We 
define a new electron operator in the rotated frame by 
a = Uc, where J7 is a 2 x 2 unitary matrix diagonalizing 
Hsd, namely, it satisfies U'^{n-<y)U = a^' . It is convenient 
to choose U ~ m cr with 



(7) 



m — I sm - cos 0, sm - sm 0, cos - ) , 



where 9 and </> are ordinary polar coordinates parametriz- 
ing the direction of n. The space and time derivative of 
the electron operator gives SU(2) gauge field 



(8) 



as dfj^c = U {d^ + i-Afj^) a, which represents the spatio- 
temporal variation of magnetization.^'* The Lagrangian 



where = (/5", j") is a four current density representing 
spin and spin-current densities ("paramagnetic" compo- 
nent) in the rotated frame, which are given by 



flta^Va- (Vat)a"a) . 



(11) 
(12) 



The spin part of the impurity potential T^mp is ex- 
pressed as Sf{c^(j°'c) ^ S°'{t){a'<a°'a), where Sf{t) = 
TZ°'^{Rj,t)Sj is the impurity spin in the rotated frame 
withi4 



n 



a/3 _ 
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(13) 



being a 3 x 3 orthogonal matrix representing the same 
rotation as U. We take a random average over impurity 
positions and for magnetic impurities we take a quenched 

average for the impurity spin direction as S*" = and 



S,j6°'^ X 



SI {a,l3 = x,y) 
SI (a, /3 = z) 



(14) 



The anisotropy axis of the impurity spin is defined in 
reference to the rotated frame. 

In the following calculation, we use the impurity- 
averaged Green's function 
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;fe + £Fa + i7crSgn(Imz) 



(15) 



where £k = h^k^/2m, epa- = sf + aM . The subscript 
a ~^,\. represents the majority and minority spins, re- 
spectively, and corresponds to ct = -1-1, — 1 in the formula 
(and to a =4,, f or —1, -|-1). The damping rate 7o- is eval- 
uated in the first Born approximation as 



7<T = 7^ = 7r(rii/„ 

2Ta 



(16) 



where v„ ~ mkpa/'^TT^h'^ is the density of states at epo 
with fcpo- = and 



n-,u 



to = 2n,u2^2 



(17) 
(18) 



with rii (ng) being the concentration of normal (magnetic) 
impurities. The first and second terms in Eq. (16) come 
from spin-conserving and spin-fiip scattering processes, 
respectively. In this paper, we assume weak impurity 
scattering, ^ epcr- 
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III. CALCULATION OF SPIN TORQUE 

Dynamics of the magnetization is described by the 
Landau-Lifshitz-Gilbert equation,^ 



h = 7o-ffcff X n + a^fi x n + r^'^, 



(19) 



Here we have introduced a four-component notation 
for the gauge potential, A™^ (-0°'", A™^), and the 

charge/ current density, = — 
with 



where 70-ffcff and ag are effective magnetic field and a 
Gilbert damping constant, respectively, in the absence of 
conduction electrons. Effects of conduction electrons are 
contained in the spin torque density. 



(22) 



jo = P = -ea)a, 

j = -eJ—{a)Va - {Va))a). (23) 
2mi 



hS 



Mn X (cr(r))r 



(20) 



which comes from Hsd- (Here, a is the lattice constant 
and S is the magnitude of d spin.) Our task is then to 
calculate (cr(r))nc — Ti-io' {r)) -ac in such noncquilibrium 
states with spin accumulation and the associated spin 
current in the presence of magnetization texture. 

To produce spin accumulation, we disturb the elec- 
tron system by applying an electromagnetic scalar po- 
tential (ff™ and a vector potential A"™ which are time- 
dependent and inhomogeneous. The perturbation is de- 
scribed by the Hamiltonian 



i?em = dr{pc 



drj^A 



cm 



(21) 



Working with Fourier components, we calculate (6-(r))r 
linear response to A 



cm 



where a = x,y takes transverse components {i.e., per- 
pendicular to the local magnetization, which is z in the 
rotated frame). The response function L'^{q, q'\ lo) is ob- 
tained from 



L^(q,q';zwA) = /'^^ dr e*"^^ (T.a"(q, r)j^(-g')) 
Jo 



l/T 



dre^"--<^(T,a"(q,T)jV(-g'))-'^, 



,'-^Y.^'^P)^^rd-^q,T)a^{-q' -p))\. 
p ) 



(25) 



r 



by the analytic continuation, — > huj + iO, where T 
is the temperature and uj\ is a bosonic Matsubara fre- 
quency. In this paper, we limit ourselves to absolute 
zero, T ~ 0. The average (• ■ ■ ) is taken in the equilib- 
rium state determined by L^i. The Fourier components 
of the spin and (paramagnetic-like) current densities are 
given by 



k 



(26) 



(27) 



with k± = k ± q/2, ~ 1 and Vi = hki/m. We extract 
an SU(2) gauge field from L'^, and write as 

L"(q, g'; *c.„) = ehL"J{q, g'; tu;„)A^{q - q'). (28) 



The coefficient L^f is expressed by the diagrams shown 



in Fig. 1 (a). 

Deferring the details of calculation to Appendices, we 
give the result here as 

1 

Yr 

(29) 



Ll^{q, q'- u + ^0) = ^{Sf (3s'^^)Kr,{q' ,u + ^0), 



where 6"^ = Sap - Sa^Spz, £xy = Syx = 1 and 



(30) 



is exactly the same as the coefRcient of the /3-terms of spin 
torque^'^"' [Eq. (1)] and spin motive force^^''^-^ [Eq. (2)]. 
In Eq. (29), Kf^ is given by 



K^,^{q,UJ + iO)^q,uK' 



K^{q, c + zO) = ^ {aDu) [ % - ^ j - ^^K^^ \ , 



(31) 

2 

(32) 
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(b) A,„= • = fif 'V)i'„+ Tsf ">t;^+ri|%A,„ + f2i.">A^s 



FIG. 1: (a) Diagrammatic representation of the coefficient 
x'^i of ^m" ^'^'i ^'^^ transverse spin polarization, 

(cr"(r))nc. The solid (double solid) lines represent electron 
lines carrying Matsubara frequency ie„ (je„ + iuj\). The 
dotted line represents scattering by impurities, either non- 
magnetic or magnetic. The filled vertex with a cross rep- 
resents the ladder-type correction to the four-current vertex 
(A^ct) as defined in (b). (b) Diagrammatic representation of 
Dyson equation for A^,,-. The dotted line with Fi (r2) denotes 
spin-nonflip (spin-flip) scattering. 



with 



(aDvY) +2ttV2{v){'7Dv) 



(33) 



Here we have defined — D^q^ — iut, {aDv) — D^v^ — 
D^v^, {aDvY) = D^v^Y^ — D^i']Y^, {v) = + v^, and 
{Yv) = Y^^v^^ + Y^u^. Note that A'f/^'s given by Eqs. (31)- 
(33) agree with the result derived in Ref. 31 and describe 
the response of spin current to the electromagnetic field, 

= e2if-(q,c. + ^0)A™^(q,c.). (34) 

Here Js('7, w) is the Fourier component of the spin-current 
density, = j'^, the right-hand side being defined by 
Eq. (12). 

Summarizing Eqs. (24), (28), (29) and (34), we have 



eh 



(8f - /3e"^) 



"° A/' 

X ^ (g - 9')A1^(g', w + *0)^7(q', w) 



^(Jf -/3£"^) 5:^^(9 -9')- Js(g', 



(35) 



From the identities, 

7e«^^f A7 = -i(nx9,n)", (36) 
7^"'^e''^A7 = -i(a,n)^ (37) 



we finally obtain the spin density as 



('T(r)). 



{nx (jVV)n-/3(j,.V)n}, (38) 



2eM 

and the spin-torque density as 

<d(^,i)-:^{a-V)n + /3nx (jVV)n}. (39) 

The first and the second terms of T'^^ir^t) represent 
the spin-transfer torque and its dissipative correction (/? 
term), respectively, induced by js given by Eq. (34). As 
shown in Ref. 31, is written as 



(40) 
(41) 



where o'er = e^Dai^a is the electrical conductivity and p„ 
is the charge density of spin-cr electrons. '^^ Therefore, js 
consists of two parts; the ordinary (local) current, ji'"-' = 
(cr-f- — <Ji)E, and the diffusive (nonlocal) current, j'i^'' = 
~{D^Vp^ — D^V pi). Equation (39) is thus a natural 
extension of the previous current-induced torque which 
contained only ji'"^ It is reciprocal to the spin motive 
force induced by magnetization dynamics"^^. These are 
the main results of this paper. 

It is noted that the spin current in Eq. (39) is within 
the two-current picture'^^"'^^ and satisfies the spin "con- 
tinuity" equation"^ ^ 



— + divj^ = 



Tal,a 



Ps 

Tal,S 



(42) 



where 



Sf ,{T 



27rF2fff/ft is the spin- flip rate for spin-cr 
electrons. Therefore, the diffusive part of the spin torque 
in Eq. (39) can be obtained by solving the ordinary dif- 
fusion equation in the two-current model obtained from 
Eqs. (41) and (42). 



IV. APPLICATION 

As an application, we consider a magnetic domain wall 
(DW) motion driven by a diffusive spin current due to 
nonlocal spin injection in a lateral spin- valve system. The 
system consists of two ferromagnetic wires, Fl and F2 
(with width wp and thickness dp), and a non-magnetic 
wire, N (with width ui^ and thickness d^). Fl and F2 
are separated by a distance L and are connected by N, 
as shown in Fig. 2. There is a single DW in F2, while Fl 
is uniformly magnetized in z direction. When an elec- 
tric current, /, is passed from Fl to the left end of N, 
as shown in Fig. 2, a spin-polarized current is injected 
from Fl into N and induces a non-equilibrium spin ac- 
cumulation in N. Then the spin current is injected into 
F2 and diffuses away from the injection point, in both 
positive and negative directions along the z-axis. The 
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motion^^ "^^ 



FIG. 2: Schematic illustration of a lateral spin- valve structure 
with two ferromagnetic leads (Fl and F2) connected by a non- 
magnetic lead (N). The bias current I flows from Fl to the 
left end of N, and produces a spin accumulation in N. The 
induced diffusive spin current flows from the Fl-N contact 
region to the N-F2 contact region and is injected into F2. The 
injected spin current diffuses in F2 in positive and negative 
z-directions away from the contact. 



spin-current density, js{z), flowing in F2 is then given by 



(43) 



where 5p,„{z) — —p„[z)/eva is the nonequilibrium part 
of the chemical potential in F2. It is obtained in Refs. 
40 and 41 as 



(44) 



where /i and 63 ^rc constants determined by the bound- 
ary condition, and Af2 is the spin diffusion length in F2 
given by'^^ 



AI2 D^D^ 



(45) 



To evaluate the effect of spin torques on the DW mo- 
tion, we assume a planar and rigid DW, a so-called one- 
dimensional model''^. The DW motion is described by 
the relevant collective coordinates, the position X of the 
DW center and the polarization angle c^o, which may de- 
pend on time. Using these variables, the DW is expressed 
as n{r,t) ~ n„{z — X{t)), where 



/ cos 00 



sm 00 ^g^j^jj^ 
ycosh(z/Aw) ' cosh(z/Aw) ' Aw 



(46) 



with Aw being the DW width. (We consider a tail-to- 
tail DW for simplicity.) In the absence of the extrinsic 
pinning, the DW motion is described by the equation of 




(47) 



sin20o+Teu, (48) 

where K± is the hard-axis anisotropy constant, and 
TVw ^ 2KA/a^ is the number of spins in the wall with 
A = (jjpdF being cross sectional area of F2 and a being 
the lattice constant. The force Fc\ and torque Td z are 
defined by42-44, 

Fei = -M J dr d^n^iz - X{t)) ■ (cr(r))„e, (49) 

Tel,, = -M fdr [nw(z - X{t)) x (cr(r))„e]_, . (50) 



By substituting the spin polarization given by Eq. (38) 
into Eqs. (49) and (50) and putting n in n^,, we obtain 



ze AwAf2 

T,u = -^^mx;0, 
Ze Af2 



where C ~ ai'^Dab'^ and 

/>oo 

J-(X;C) =2^ / dze-l"+^l/« 



sinh z 
cosh'^ z ' 



(51) 
(52) 

(53) 



with X = X/Aw and ^ = Af2/Aw As pointed out in 
Ref. 41, a large spin-current injection from N into F2 is 
possible in a device with a tunnel junction at the injection 
part and a metallic contact with a strong spin absorber 
at the detection part of the nonlocal spin valve structure. 
In such a case, we have 



C: 



A3 



(54) 



where Aj = ojnWf is the contact area of the junction 
and Pi is the spin polarization of Fl and An is the spin 
diffusion length in N. Thus we obtain 



ze Aw 



where 



el, 2: 



A 



2e 



(55) 
(56) 

(57) 



is the effective spin current due to spin diffusion. Fig- 
ure 3 shows J^{X;^) as a function of X for several val- 
ues of ^. Two (positive and negative) peaks are seen at 
around \X\ ~ 1 for each ^, which corresponds to the dis- 
tance between the wall center and the injection point. 
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FIG. 3: Graph of T{X;^) as a function of X = X/X„ for 
several choices of ^ = Af2/Aw 



This is due to properties of that diffusive spin current 
flows in gradient direction of the spin accumulation. In 
a one-dimensional wire, the diffusive spin current flows 
in positive and negative z-directions and drives the wall 
in mutually opposite direction separating the injection 
point. Note that its direction can be reversed by revers- 
ing the current direction without reversing the magne- 
tization of Fl, or by reversing the magnetization of Fl 
without reversing the current direction. We also see that 
the force and torque acting on the DW strongly depends 
on the values of ^. For smaller ^, the peak value of J-' is 
also smaller, because of the shorter spin-diffusion length. 
These facts show that an efflcient angular momentum 
transfer can be realized if we could generate the spin ac- 



J 



cumulation at the around edge of the DW and we use 
materials with large value of ^, where the spin diffusion 
length Af2 is much longer than the wall width. 



V. CONCLUSION 

We have presented a microscopic theory of current- 
induced spin torque with emphasis on the role of diffusive 
spin current. The obtained torque consists of reactive 
and dissipative terms, and provides a natural extension of 
the well-established one due to 'local' spin current which 
accompanies ordinary (Ohmic) electric current. We ap- 
plied the result to a domain wall in a lateral spin-valve 
system with nonlocal spin injection. The effects on the 
domain wall (force and torque) depend strongly on the 
ratio of the spin diffusion length to the wall width, and 
increase when the spin diffusion length is much longer 
than the wall width. 

An interesting extension of the present theory would be 
to the study of modulation of the magnetization damping 
due to spin accumulation, as observed experimentally.'*^ 
This will be reported elsewhere. 
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Appendix A: Calculation of L 



a/3 



In this Appendix, we evaluate the response function L"f in Eq. (28). It is expressed as (see Fig. 1) 

L"J{q, g'; ^u;,) = tJ2J2{''^ ^"^+^^-,72 ^t'^"GU,'/2 KGu-^' 12 ] 

+ tr[S"G+^,V2^MGfe-,72-r^"G,_,+,v2] + ^tr[I]"G+^^/,a^Gfe_,/^ (Al) 

where vf = {h/m) {ki ± qi/2), and we have adopted a 2 x 2 matrix notation, (G)o-o-' = GcSaa' ■ We have included a 
ladder correction to the four-current vertex via A^ ^ + A^ (A^ being deflned in Fig. 1(b)), and a first-order 
correction to the transverse-spin vertex, tr" (a = x,y), via 



(A2) 



where Fq = riiu^ — naU^S"^. In the following calculation, we consider the adiabatic limit and put q = in the Green's 
functions. Namely, we consider L^f (q', zwa) = i^f (0, q'; iw^), or using q for q'. 



fit 



^ E E A^Gk- ] + tr[ I]"G++ A,,Gfc_ v,<j^Gk+ ] + ^ tr[ I]"G++ a^Gk- ] | , 

n k ^ ' 



(A3) 
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where fe± = fc ± q/2. Note that, here and hereafter, q actuahy means q' of We retain terms up to 0{uj, q^), 

and to the next leading order with respect to the electron damping ja- to sec the spin-relaxation effects. Equation 
(A3) is written as 

L''J{q,ioJx) - (5f +ia£"^) ^^,,,{q;ien + iLOx,ien), (A4) 

n (T 

with 

k ^ m J 

(A5) 

where 

r. = fo^G+^Gfc.. (A6) 
Performing the analytic continuation, iuj\ — > hio + iO and retaining terms up to the first order in w, we obtain 



L;f(9,.. + zO)-i^f(9,0) 



(A7) 



where e± ~ e ± uj/2 and /(e) is the Fermi distribution function. The superscripts (1) and (2) indicate the analytic 
continuations, G(ie„ + iLox)G{ien) -> G^G^ and G^G^, respectively. The ^p''^^^ and (/^^^'^ are given by 

^i.t(g;0,o) = (i + rw) |i3,_(q) + s^,,,(-q) + i|-rw|, (as) 

'pS-(9;^,0)= (l + r(2)) |G^.,.(q) + G;,,^(-q) + A(ing,^)Co.,.(g)+A(^j(q,a.)GoV(-9) + ;^ (A9) 
where 

B^^.A^l) = E «m^'-GL,^G«+,,GL,,, (AlO) 
fc 

G^.,.(g) = ^ «^tvG^,^G«+,,Gt,,. (All) 

Since the terms of our interest are 0(7^) or lower order in 7^, the four-current vertex correction A^ is relevant only 
in (p(^) and has been dropped in ip (1). The A(,l^(q , uj) is evaluated in the diffusion approximation by retaining q and 
a; in a usual way, with a result 



31 



where Y^ = Do-g^ - iijJ, {v) ^ vi^ + i^i, {Yv) = Y^v^ + Y^v^ and {Dv) = D^^v^ + D^v^, with 1?^ = 
(vpa = hk-pa/Tn) being a diffusion constant. 

In B^^ u and G^v_cr, Q is retained up to the first order whereas w is set to zero, lu ^ 0. These are calculated in 
Appendix B, and the results lead to Re(/3^^j'g. = 0(70-) and 



^!fi(g; ^' 0) = -^(1 + ^'^Z^)^ (nil + nloU^^)) + ^(i + *^/5)a (n^^^, + u^.f^ 



lira 



+ (A13) 
where 

n|,ll. = E ^^M^.G«+.,Gt^,. (AM) 
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The last term in the braces of Eq. (A7) is calculated as 



(A15) 



which cancels the terms in the second line of Eq. (A13). Therefore, we obtain 

L^f (g, . + .0) - L^f (g, 0) = ^(^f - Z?^"'^)^ E ^ ("£1 + 



or, in the original notation in Eq. (Al), 



L^f + *0) ^ L;^(O,q';0) + -(<5f - + zO), 



(A16) 



(A17) 



where 



i^-(g, c + .0) = ^ ^ a (n(^)^ + n(^)^A(S) (A18) 



is nothing but the electromagnetic response function of 
spin current js^i (see Eq. (34)) obtained in Ref. 31. The 
first term of Eq. (A17) gives a spin torque due to dia- 
magnetic spin current, which will be reported elsewhere. 



Appendix B: The fc-integrals 

The k integrals of the product of Green's functions 
are given here. We retain low-order contributions with 
respect to 7^ and q. 

Equation (A14) is calculated as'^^ 



The Cij,a{q) in Eq. (All) is calculated as 

k 



^ , ■ Tt(2) , A 



(B4) 



where 274- = 7-t- + 74., and ps ~ — is the spin 
density with Uc — l£-paV„ ji. In the last line, we have 
used Eq. (A14). 

The Cio,cr(g) in Eq. (AlO) is calculated as 



(2) 
00,0- 

(2) , 
iO,cr - 

(2) . 



2t: D^VcrSij, 



(Bl) 
(B2) 
(B3) 



2{aM - ^7+) 



'^iGk+.a{Gk-,a ^ G'^_^) 



k 



2AI V mJ '^-^ 



(B5) 



where = Vp^Ta-/3 with upo- ~ hk-pa-/m. 



Finally, Eq. (A9) is calculated, with E^^' = -iTrarov+/2M + 0(7), as 

(i + r(,2))(c^,,^ + A(,i)Co.,.) 



2M \ 2M " +A M 



(2) 



4mA'f • 



where 



lira 



(B6) 
(B7) 
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This leads to Eq. (A13). 
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